Towards a General Theory of 

Simultaneous Diophantine Approximation of 

Formal Power Series: 

Linear Complexity of Multisequences 

)0 ■ Michael Vielhaber 1 and Monica del Pilar Canales Chacon 1 

o 
o 

Institute* de Matematicas, Universidad Austral de Chile, 
Casilla 567, Valdivia, Chile {vielhaber ,monicadelpilar}@gmail . com 

Abstract 

We model the development of the linear complexity of multise- 
quences by a stochastic infinite state machine, the Battery-Discharge- 
Model, BDM. The states s G S of the BDM have asymptotic probabil- 
c*2 . ities or mass Hoo(s) = V(q, M)~ l -q~ K ( s \ where K(s) G No is the class 

of the state s, and V(q,M) = £^ eNo PM(K)q~ K = Uti <?7(<f " 1) 
is the generating function of the number of partitions into at most M 

■ parts. We have (for each timestep modulo M + 1) just Pm{K) states 
O ! of class K. 

We obtain a closed formula for the asymptotic probability for the 

■ linear complexity deviation d(n) := L(n) — \n ■ Mj{M + 1)] with 
O ' 

o 

Zj . The precise formula is given in the text. It has been verified numeri- 



X 
S3 



7 (d) = e (V M(M+1) ) , VM G N, yd £ z. 



cally for M = 1, . . . , 8, and is conjectured to hold for all MgN. 

From the asymptotic growth (proven for all M G N), we infer the 
Law of the Logarithm for the linear complexity deviation, 

. , d a (n) 1 d a (n) 

— iim mi = —— — — = iim sup , 

n^oo logn (M+l)logg n ->oo log n 

which immediately yields a ^ - — > with measure one, VM G N, 

n M + 1 

a result recently shown already by Niederreiter and Wang. 

Keywords: Linear complexity, linear complexity deviation, multi- 
sequence, Battery Discharge Model, isometry. 
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1. Linear Complexity of Multisequences 

The linear complexity of a finite string a G F™, L a (n), is the least length of 
an LFSR (Linear Feedback Shift Register), which produces a\, . . . , a n start- 
ing with an initial content a±, . . . , o,L a (n)- If all symbols are zero, we set 
£(o,...,o) (n) = 0. Also, we put L a (0) = for all a. 

An alternative and equivalent definition defines L a (n) as the length of the 
shortest recurrence within the Oj, i.e. 



i-i 



L a {n) := min 3a u . . . , ct;_i G F ? , VI < fc < n - I : a fe+ ^ = V] a* • a k+i 



Given an infinite sequence a G F^°, we define L a (n) as before, taking 
into account only the finite prefix ai,...,a n . The sequence (L a (n)) neN is 
called the linear complexity profile of a. The diophantine approximation of 
the generating function G(a) := Y^=i a n x ~ n G F g [[a; _1 ]] by a polynomial 
function with precision at least k that is 

G(a) = ^\ + o(x- k ), 
v(x) 

requires a polynomial v(x) of degree at least L a (k), and this length is also 
sufficient, since v(x) may be chosen as the feedback polynomial of the LFSR 
producing a±, . . . , a^. 

Turning to multisequences (a n ,m)neN,i<m<M G (F^) , we ask for simul- 
taneously approximating all M formal power series 

OO , s 

F.tfx- 1 ]] 9 G m (a) :=J2 a n,mX- n = + o(x- fc ), 1 < m < M 

vyx) 

n=l v ' 

with the same denominator polynomial v(x), equivalently we search a single 
LFSR which produces all M sequences with suitable initial contents. The 
linear complexity profile of a now is defined by (L a (n, m))n€N ,i<m<M (symbol 
by symbol), and we set L a (n) := L a (n, M), when considering only complete 
columns of all M sequences at the same place n, with profile (L a (n)) ne ^ . 
The goal of this paper is to characterize the behaviour of L as a 
probability distribution over all multisequences from (F^)°°. 
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2. Continued Fraction Expansion: 
Diophantine Approximation of Multisequences 

The task of determining the linear complexity profile of one multisequence 
from (F^) 00 has been resolved by Dai and Feng [2]. Their mSCFA (multi- 
Strict Continued Fraction Algorithm) computes a sequence 



(n,m)e{l,...,M}xN 

of best simultaneous approximations to ((?*•). The order of timesteps is 
(n, m) = (0, M), (1, 1), (1, 2), ... , (1, M), (2, 1), (2, 2), . . . with 

Gm(a) = ^2 

Q"m,t ' % 

We will denote the degree of v (jl ' m \x) by deg(n, m) G No, thus the linear 
complexity profile is (deg(ra, M)) neNo = {L {Gmil < m < M) (n)) nmQ . 

The mSCFA uses M auxiliary degrees W\, . . . ,wm G No- The update of 
these values (and deg) depends on a so-called "discrepancy" 5(n,m) G ¥ q . 
S(n, m) is zero if the current approximation predicts correctly the value a nim , 
and S(n, m) is nonzero otherwise. 

Furthermore, the polynomials u m (x) and v(x) are updated, crucial for the 
mSCFA, but of no importance for our concern, and we omit the respective 
part of the mSCFA in the program listing: 

Algorithm 1. mSCFA 

deg := 0; w m := 0, 1 < m < M 
FOR n := 1,2, ... 

FOR m := 1, . . . ,M 

compute S(n, m) / / discrepancy 

IF 5{n,m) = 0: {} // do nothing, 0, Thm. 2, Case 2a] 
IF 5(n, m) ^ AND n - deg -w m < : {} // Thm. 2, Case 2c] 
IF sin, m) ^ AND n - deg -w m > 0: // Thm. 2, Case 2b] 
deg_copy := deg 
deg :=n-w m 
w m :=n- deg_copy 
ENDF0R 
ENDF0R 




(n.m) / \ 



,(n,m)i 



X 



+ o(x~ n ),ne N ,V 1 < m < M. 
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3. The Battery-Discharge-Model 



This section introduces the Battery-Discharge-Model (BDM), a stochastic 
infinite state machine or Markov chain, which will serve as a container to 
memorize the behaviour of deg in the mSCFA for all inputs a e 
simultaneously. 

Since the linear complexity grows approximately like L a (n) = deg(n, M) ~ 
\n ■ jjfW] (exactly, if we had always S(n,m) ^ 0), and the auxiliary degrees 



w m of the mSCFA grow like w n 
average behaviour as follows: 



l_ IT+iJ ' we ex ^ rac t the deviation from this 



Definition. The linear complexity deviation or degree deviation is 

M 



d := d a (n) := deg- 



n 



E Z, 



(1) 



M+l 

which we call the " drain" value, and the deviation of the auxiliary degrees is 

1 



b m : = 



n 



M+l 

which we call the " battery charges" . 



w m G Z, 1 < m < M, 



(2) 



The BDM will assemble all necessary information about the development 
of d a (n) with n — > oo and a running through all of (Ff)°°. 

We establish the behaviour of d and b m in two steps. First we treat the 
change of d, b m when increasing n to n + 1 in the mSCFA (keeping deg, w m 
fixed for the moment): 



deg- 

and 
(n + 1 



(n + 1) 



1 



M 



M + l 



M+l 



W m = 



deg - [■ 
deg - [■ 



M 
M+l 



- 1, n^MmodM+1, 
n = MmodM +1, 



n ^ MmodM + 1, 



l n ■ mTi J ~ w ™ + X > n = MmodM +1. 



(3) 



(4) 



Hence, by (3) we have to decrease d in all steps (we call this an action 
gL), except when n = M — > n = mod(M + 1), and only here we increase 
all M battery values 6 m , by (4) (action b+). 
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With d(M, 0) = b m (M, 0) := 0, Vm, we obtain the invariant 



d(n, m) + 



'b k (n,m) + nmod(M + 1) = 0, Vn G N , 1 < m < M (5) 



for the initial timestep (n, m) = (M, 0). Also, by (3) and (4), the actions oL 
(increase n, decrease d) and 6 + (decrease n mod (M+ 1) by M, increase M 
batteries by 1 each) do not change the invariant. 

Now, for n fixed, the M steps of the inner loop of the mSCFA change w m 
and deg only in the case of S(n, m) ^ and n — deg — w m > that is 



n — deg — w m > 



n 



d + 



M 



n ■ 



M+l 



n 



M + l 



-6m > 



b m > d. In the case 5^0 and b m > d, the new values are (see mSCFA) 



deg = n — w m and 
and thus in terms of the BDM variables: 



d + ( = 6) (n - w m ) - 



w. 



n 



deg 



(6) 



" n ■ M ' 


(2) 


n 


+ b m - 


n 


M+l 




_M+1_ 


_M+1_ 



and 



(2;6) 



n 



M+l 



- (n - deg] 



(i) 



n ■ M 
M + l 



+ d + 



n-M 
M + l 



d. 



an interchange of the values d and b m . We say in this case that "battery b m 
discharges the excess charge into the drain" , and call this behaviour an action 
"D" of battery b m , corresponding to case 2b of Thm. 2}. A discharge does 
not affect the invariant (5), which is thus valid for all timesteps (n, m). 

The remaining cases are b m > d, but 5 = 0, an inhibition of b m , action 
"I" (case 2a of Thm. 2]), and two actions of do nothing, "N = " and "N<", 
distinguishing between b m = d and b m < d (case 2c and part of CcLSG 2cl). 

Since we do not actually compute the discrepancy 5 (in fact, we do not 
even have a sequence a), we have to model the distinction between 5 = and 
5^0 probabilistically. 

Proposition 2. In any given position (n,m),n e N, 1 < m < M of the 

formal power series, exactly one choice for the next symbol a n ^ m will yield a 
discrepancy 5 = 0, all other q — 1 symbols from ¥ q result in some 5^0. 
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Proof. The current approximation Um (x) / v^ n ' m \x) determines exactly 
one approximating coefficient sequence for the m— th formal power series 
G m . The (only) corresponding symbol belongs to 5 = 0. □ 

In fact, for every position (n, m), each discrepancy value S G ¥ g occurs 
exactly once for some a„ im G ¥ q , in other words (see PO^II] for M = 1): 

Fact The mSCFA induces an isometry on (F^) 00 . 

Hence, we can model 5 = as occurring with probability 1/q, and 5^0 
as having probability (q — l)/q. 

To keep track of the variables d, b m , we define the following state set for 
the BDM: 

Definition. The augmented state set is 

S := {s = (bi, . . . , b M , d; T, t) | b m G Z, 1 < m < M; d G Z; 

M 

T G Z; 1 < t < M + 1; d + T + ^ 6 m = 0}, 

m=l 

where the last condition is the invariant (5). For the BDM, we only use the 
timesteps < T < M, and the BDM thus has the state set 

S := {s G S | < T < M} 

with initial state s := (0, . . . , 0; 0, M + 1). 

To facilitate notation, we also define S(T ,to,do) = {s G S \ T(s) = 
T ,t(s) =t ,d(s) =d }, and similarly S(T , t ), S(T ), S(T , t , d ). 

A state stores the values of the batteries and the drain in b\, . . . , bM, d, 
the value T corresponds to the time modulo M + 1 that is T = n mod M + 1, 
and the "ministeps" t = 1, . . . , M correspond to the update of battery b m 
between t = m and t = m + 1, while t — M + 1 — > 1 corresponds to the 
updates d^,b + . 

The allowed transitions a (action) from a state s = (bi, . . . , 6m, d; T, t) 
are a = gL or 6 + for £ = M + 1 , and otherwise depend on the relative size of 
b t and d (a G {-D, /, N=, N K }). We have s — > s + with the following actions, 
conditions, nextstates s + , and probabilities: 
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a 
D 
I 

N = 



d. 

&4 



condition 

b t > d,t< M 
b t > d,t<M 
b t = d,t<M 
b t <d,t<M 
I < M.I .1/1 
T = M, t = M + 1 



(.01, 
(h, 
(h, 
(h, 



n, 



prob. 
q - 1/q 
1/q 
1 
1 
1 
1 



, bt-i, d, bt+i, . . . , 6m , W, T,t + 1) 
,b M ,d; T,t+1) 
,b M ,d; T,t+1) 
,b M ,d; T,t+1) 
,b M ,d- 1;T + 1,1) 
(6i + l,62 + l,...,6M + l,d;0,l) 
Whenever b t > d, both D and I may occur, leading to two feasible tran- 
sitions from a given state s, whose probabilities sum up to 1. 

Recall that from (T, M + 1) to (T + 1, 1), the drain d is decremented 
according to (3) for T < M, action "d_", and from (M,M + 1) to (0, 1), the 
batteries b m are incremented according to (4), action "b + ". 

Definition. The state transition matrix T of the BDM is an infinite stochas- 
tic matrix indexed by s, s' G S, and where 



T(s,s') 




(q-l)/q, s 



D 



s , s — > s , s 



s , or s — >■ s , 



otherwise. 



Every row either includes an "I" and a "D" , or else one of "N = 



"d_", or "b + ". Reading the feasible transitions backwards, one obtains that 
a state with b t < d (at (T, t + 1)) is reached either by a discharge, or by a 
"N<", hence the corresponding column of s' sums up to 2— ^ + 1. Likewise, if 
s' has b t > d, this may only be the result of an inhibition, hence column sum 
1/q. The cases "N=", "d_", and "b + " all are by themselves the only nonzero 
entry within a column, which has thus sum 1. 

In terms of d, b m , we have the following equivalent probabilistic formula- 
tion of the mSCFA (timestep t — M + 1 comes after the FOR m = t loop): 

Algorithm 3. BatteryDischargeModel 

d := 0;b m := 0, 1 < m < M 
d := d — 1 1/ gL 
FOR n := 1,2, ... 

FOR m := 1, . . . ,M 
IF b m > d: 

WITH prob. (q - l)/q: 

swap(6 m ,c?) // action D 



7 



WITH prob. 1/q: 
{} I) action / 

ELSE 

{} // action iV=,iV< 
END IF 
ENDFOR 

IF n ^ MmodM + 1 : 

d:=d-l II d- 
ELSE 

b m := b m + 1, 1 < m < M If b. 
END IF 
ENDFOR 



4. Classes of BDM States 

The Markov chain BDM will turn out to be strongly concentrated on few 
states. We define a family of measures \x T on S, indexed by r G No. We start 
for r = with all mass concentrated on the initial state s : 



1, s = s = (0,...,0;0,M+ 1) 
0, s ^ s 



For successive timesteps r, we then put /x r+ i(s') = /i r (s), if s s', 
s —> s', or s -^U s'. Also, /i r+ i(s') = -// T (s), if s — s'. Finally /i r+1 (s') = 

^Vr(si) + A*t(s2), if si — s' and s 2 ^> s'. 

Put in other words, (/i r +i(s)) se5 = 7" • (A i r(s)) se5 - 

Be aware that from ji T to /i T +i, we only deal with one input symbol (or 
gL, hence the distribution after reading all M inputs of column n is in 

fact /!(Af+l). n (s). 

Definition. We will use repeatedly the "timesteps" (T,t) G {0,...,M} x 
{1, . . . , M + 1} of the BDM, comparing them with linear time r € N . We 
define: 

(T, t) = r :<^> (T - 1) • (M + 1) + t = t mod (M + l) 2 

When dealing with the m-th symbol in column n, the r-th input symbol, 
we are in a state with T(s) = n mod (M + 1), t — m, and (T, t) = r. 
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Proposition 4. For every r G No, 

(s) = { !> ( T o,*o) = r 

Proof. By induction on (T, £): Initially (r = 0,T = 0,t = M + l),all mass is 
on So- Also, every transition goes from states with (T, t) = r to states with 
(T', £') = r + 1, carrying over the mass to the new S(T', t'). □ 

Definition. Denote the number of sequence prefixes in (F^)" with linear 
complexity deviation (IgZas N(n, d; q). 

Since the BDM has been derived from the behaviour of the mSCFA, we 
obtain 

Theorem 5. Assume that exactly N of the q M ' n sequence prefixes of length 
n lead to a certain configuration ( deg, w±, wm) of the mSCFA, and let 
bi, . . . , 6m, d, T be derived from n, deg, w m according to (1), (2), then 

N 

H( M +iyn(bi, ■ ■ ■ , b M , d; T, M + 1) = -p— 

with (T, M + 1) = (M + 1) • n that is T = n mod (M + 1). 

Proof. The theorem is true for n — 0, (T,t) = (0, M + 1) with N — 1, 
deg = d = b m = w m = 0, Vm, starting with the (only) prefix e, the empty 
string. 

From then on, by the construction of the BDM, for t(s) < M a transition 
s — > s* takes place with probability ^, with 6 from {1, g — 1, g}, if and only if 
the mSCFA goes to the state corresponding to s* for b out of the q possible 
next symbols a n)Tn G W q , or, for t(s) — M + 1, t(s*) = 1, corresponding to 
actions d-,b + , with probability one, while the mSCFA increases n. □ 

From this theorem now follows as a corollary the description of N(n, d; q) 
by the mass distibution on the BDM states (more on finite n in Section 9): 

Theorem 6. For n G N with (T , M + 1) = [M + 1) • n, for d G Z, 

iV(n, d; q) = q M ' n X /i(M+i)-n(s) 
seS(T ,M+l,d) 



9 



Definition. For a given state s G S, we define its asymptotic measure as 
//oo(s) := lim sup /i„(s) = lim /i„(s). 

n=(T( S ),t(«)) 

We have ^ se5 A*oo(s) = (M + l) 2 , since each S(T,t) sums up to 1. 

We will see that all states satisfy pLoo(s) = ^oo(s )-q~ K ^ for some K(s) G 
No. We call this value K(s) the class of state s and define it algorithmically, 
generalizing to s G S: 

Definition. The class of a state s = (pi, . . . , 6m, d; T,t) E S is 

M+l 

tf(s) = -tt s + M ■ T + 2 • ^ 6 m • (M + 1 - m), 

m=l 

where 7r s is minimum number of transpositions between neighbours necessary 
to sort (bi, . . . , bt-i, d, b t , . . . , 6m) into decreasing order as (b\, . . . , &m+i), h > 
h+i, 1 < % < M). Observe that the place of d in the initial sequence depends 
on t. 

Example. The state s = (-5,4, -4,2; 1,2) with M = 3, d = 2, T = 1 and 
t = 2 requires the sorting of (—5,2,4,-4) into (4,2,-4,-5), using ir s = 4 
transpositions, and thus K(s) = -4+3-l+2(4-3+2-2+(-4)-l + (-5)-0) = 23. 

This static way of determining K(s) is compatible with the following 
dynamic consideration of transitions. First we need a technical lemma: 

Lemma 7. For all s = (bi, . . . , 6m, d; T, t) G S, we have 

K(s) = K(b ± , ...,b M ,d; T, t) = K(b x + 1, . . . , b m + 1, d + 1; T - M - 1, t). 

Proof. Let s = (b u . . . ,b M ,d;T,t) and s' = (b[, . . . ,b' m ,d';T',t) := (6i + 
1, . . . , b m + 1, d + 1; T — M — 1, £). Then 7r s = 7r^, since the relative order 
within s and s' are the same, also — bi + 1,1 < i < M + 1. Using 

2-^l-(M + l-m)=2( ^ )=M-(M+1), 

m=l ^ ' 
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we thus have 



M+l 

K{s) = -n s + M-T + 2-J2~ b m-{M+l-m) 

m=l 

M+l 

= -n' s + MT- M(M + 1) + 2 • ^(6 m + 1) • (M + 1 - m) 

m=l 

M+l 

= -n' s + M -T' + 2- ^V m -(M + l-m) = K(s'). □ 

m=l 

We now obtain the change in class by counting actions / and N K : 
Theorem 8. 

(i) For every feasible transition s — > s' between states s,s' G S with 
a G {D,I,N=,N < ,d-,b + }, we have 

( 1, a = I 

K(s') = K(s) + 1 -1, a = iV< 

{ 0, a G {-D, N=, d-, b + } 

(ii) Let s a — " k s be some path from the initial state s to s. 

Let #1 = #{1 < i < k \cti = 1} and #iV< = #{1 < i < k \oi = N<}. Then 

K(J)-K( 8 )=#I-#N < . 

Proof, (i) We first deal with transitions from D, I, N = , N < that is t ^ M + 1. 
Since the values of the multiset {bi, . . . , &m, d} only get swapped (in the case 
of a discharge D) , the sum 2 J2 m =i • (M + 1 — m) as well as the term M ■ T 
stay the same. 

It suffices thus to compare tt s with tt s i. Let (b±, . . . , b t -i, d,b t , . . . , b^/) be 
the sequence to be ordered for s and similarly for s', after the discharge, 
(b[, . . . , b' t , d', b' t+1 , . . . , b' M ). Since b m = b' rn for m ^ t, they need the same 
number of transpositions, and we may in fact restrict our comparison to the 
sorting of (d, b t ) for s and (b' t , d) for s' . 

Case a = D: We had b t > d (to be able to apply D), d! := b t ,b' t := d 
and thus d' > b' t . Both before and after the discharge, one transposition is 
necessary and thus n s = 7iy, K(s) = K(s'). 

Case a — I: Again b t > d, b' t := b t ,d' := d. We sort (d,b t ) with one 
transposition, (b' t , d') is already sorted. Hence tt s > — tt s — 1, K(s') = K(s) + 1. 

Case a = N = : Now, b t = d = b' t = d' and so n s r = n s , K(s') = K(s). 
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Case a = iV<: Here b t < d and b' t < d' . (d, bt) is already in order, (b' t , d 1 ) 
requires a transposition and so ir s i — 7r s + 1, K(s') = K(s) — 1. 

Finally, for oL we have s' = (bi, ... , &m, d— 1; T + l, 1) from t = M + 1 for 
s. We have to sort (61, . . . , &m, cf) for s (with t = M+l) and (d—1, £>i, ... , 6m) 
for s' (with t = 1). Sorting first the part 61, ... , bu with equal effort in both 
s and s', 7f s = W s /, we then introduce d from the right, respectively d—1 
from the left to the same place: (61, ... , 6fc_i, d or d — b^+i, . . . , 6m), where 
fye-l > > d — 1 and bk+i < d — 1 < d. The total number of transpositions 
is then n s = W s + (M + 1 — k) and 7r s / = W s / + (k — 1). 

The class now is 

M+l 

K(s) = -tt(s) + M • T + 2 ^ 6 m (M + 1 - m) + 2rf(M + 1 - jfe) 
= -tt(s) + (M - (jfe - 1)) - (Jfe - 1) + M ■ (T + 1) 

A/+1 

+2 b m (M + l-m) + 2(d-l)(M-(k-l)) = K(s'). 

m=l,m^k 

The case a = b + is equivalent to gL followd by incrementing all the b m 
and d, hence follows from the case a = gL and Lemma 7. 

(22) This follows by applying (i) to the /fe transitions leading to s, starting 
in so with K(sq) = 0. □ 

We will now show that the limit mass distribution fi^ follows in fact (up 
to a constant) from the state classes as /ioo(s) = Co • q~ K( > s \ First, we state 
a theorem by Rosenblatt (an infinite matrix version of Perron-Frobenius): 

Theorem 9. (Rosenblatt, Let T be a Markov chain, finite or infinite. 
"If the chain is irreducible and nonperiodic, there is an invariant instanta- 
neous distribution if and only if the states are persistent, in which case the 
distribution is unique and given by [8, p. 56], where Uj = lim^oop^, 

and p£j is the probability to return to state j after n steps. 

Proof. See jHl p. 56]. □ 

Here T certainly is periodic, with period (M+l) 2 . The (M+l) 2 -th power 
of T has the property that transitions occur only within the sets S(T,t), so 
it can be ordered into a block diagonal matrix. We use only the block with 
(T,t) = (0, M +1), including s , as f := T( M+1 ) 2 \ seS{ o,M + i)- 

T and thus T is irreducible, since we get from sq to every state and back 
by the following theorem: 
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Theorem 10. (i) For every state s G S, there is exactly one sequence of tran- 
sitions a — aii • • • ex k with s =Q s and a(s) G {D, I, N=, gL, b + }* (avoiding 
actions of the type iV<) 7 which touches the state Sq only initially, 
(ii) Also, there is exactly one path from s to s avoiding actions of type I, 
which touches the state s only finally. 

Proof, (i) Unicity: There is at most one such transition: When going back- 
wards from s to So , running through the batteries in reverse order M, M — 
1, . . . , 1 for each transition, we have: 

For b m < d this results either from a discharge D, or a do nothing N < . 
Since iV< is not allowed in a s , put D. 

For b m = d this results from a do nothing N=. 

For b m > d, only an inhibition I is possible. 

Existence: There is an infinite chain of predecessors, all of class less than 
or equal to K(s). Since for each K, there are only finitely many states with 
this class, in particular, there is some state s* with (T,t) = (0, M + 1), 
which is reached repeatedly. If this state is s* = s , we are done. If not, 
s* = (bi, . . . , 6m, d; 0, M + 1) has mx := max{6i, . . . , bu-, d} > 1 and mn := 
min{6i, . . . , bu, d, } < — 1 (by the invariant (5), with T being zero). 

However, a cycle s* ^ ' - — > + ^ s* without / or ^< is impossible: Either 
d = mx at (T, t) = (0, M + 1), or else some battery b t = mx has to discharge 
(I prohibited). At (T,t) = (1,M + 1), we have d = mx > 1 in any case, 
thus at (T, t) = (2, 1), we get d > 0. Now, since mn < —1 is the value of one 
of the batteries, say bt*, at time (2, t*) we have bt* < d and thus iV< is the 
only possible action. So, no return to s* avoiding I and N K (having reached 
K — 0, there is no further decrement) is possible, unless s* = s . Since 
the only cycle to avoid passes repeatedly through s , a(s) is well-defined by 
excluding this case. 

(ii) To get back, just choose D, whenever b t > d. In this way, the class 
can never increase, and thus eventually, we must hit a cycle. But we have 
already seen that the only cycle avoiding both I and iV< passes through the 
states with class 0, including s . □ 

Theorem 11. For any two states s, s' G S, 

^oo( S ) _ K(s')-K(s) 

Proof. Let a mass distribution /i(s) := q^ K ^ be given. We show that /i 
is invariant under the transition matrix of the BDM, i.e. (/i(s)) se s is an 
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eigenvector of eigenvalue one, and unique with this property up to a constant 
factor. We consider all states leading to a fixed state s. We have three cases: 

1. b t < d after the action, coming from S\ s or s 2 — ^ s, and thus 
/-too(s) = ^-/ioo(si) + ^00(82). Since N < decrements the class (but D not), 
we have ^-p Q ■ q~ K ^) + Po ■ q ~ K ^) = ^- Vo ■ q ~ K ( s ) + Po . q -(K(s)+i) = 

(^ + l)-Po-q~ K{s) =Po-q~ K(s) - 

2. bt — d after the action, which must be a do nothing, a = N=, and thus 
K(s) = K(s'), /ioo(s) = /ioo(s')- 

3. b t > d afterwards (and before), from an inhibition, a — I which 
increments the class, hence q~ K ^ ■ ± = q-( K ( s )+ 1 ) 

This shows consistency of fj,(s) = c ■ q~ K ^ with the behaviour of the 
BDM, or stated otherwise: (/z(s)) se s = (q~ K ^) s es is an eigenvector of the 
infinite state transition matrix of the BDM. Furthermore its eigenvalue 1 is 
the largest eigenvalue of T, since T is stochastic. 

Now, T inherits the eigenvector /x, restricted to states from S(0, M + 1), 
with eigenvalue l( M+1 ) = 1. This matrix is aperiodic and irreducible by 
Theorem 10, and by Theorem 9 (Rosenblatt), fi is already the only such 
eigenvector up to a constant factor, and it remains to normalize it. 

Returning from T to T, we obtain the statement, since //(s) = q~ K ^ for 
all s G 5(0, M + 1) forces all other states in 5 also into this eigenvector. □ 

5. Antisymmetry 

In this section, we consider only the configurations with t — M + 1, at the 
end of a complete column from the input a. 

Proposition 12. For M e N, T e Z, d e Z, k e N , and 2 < q e N, Zei 

A = {s E ~S{T, M + 1, d) I X(s) = A;} and 

A = {se S(M - T, M + 1, -d) I K(s) = k}. Then \A\ = \A\. 

Proof. We show that states s := (b u ... , 6 M , — T — X; T, M + 1) and 
s := (-6 M - 1, ■ ■ ■ , -61 - 1, T + X; M - T, M + 1), where X := £™=i 6 m , 
satisfy: 

(i) T(s) = T and T(s) = M-T, (ii) d(s) = -d(s), and (Hi) K(s) = K(s). 
Then s E A sG i, and we have a bijection between A and A, hence 
\A\ = \A\. 

(i) and (ii) are obvious by inspection. To show (Hi), we first sort bi, . . . , bu 
of s into decreasing order as b\ > b 2 > • • ■ > b M by n' s permutations of neigh- 
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bours. Then —6m — 1, • • • , ~bi — 1 of s can be sorted into decreasing order 
by tt' s permutations at the same places into — &m — 1) • • • , ~b\ — 1. 

We now introduce d = — T — X and — d = T + X, resp., into the ordered 

6's as bi > . . . b k > —T — X > b k+ \ > . . . b M and — b M — 1 > b k+ i — 1 > 

T + X > —b k — 1 > . . . bi (observe the > inequality in both cases to the 
right of ±(T + X)). We have a total of 7r s = tt' s + M - k and ttj = < + fc, 



M 



resp., permutations, thus K(s) = -n s -(M-k)+M-T + 2 J^(M + 1 - 

i=i 

M 

i)6i - 2 bi + 2(M - fc)(-T - X) and K(s) — —n s — k + M ■ (M — T) + 

i=k+l 

M k 

2 J^(M + 1 - i)(-b M +i-i - 1) - 2 J^(-6i - 1) + 2A;(T + X), where the first 

i=l i=l 

sum treats the fej's in their place &e/ore introducing ±(T + X), the second 
sum adjusts the foj's, which are shifted while introducing d, by 2, and the last 
term belongs to the drain ±(T + X). The difference is then 

M 

K(s) - K(s) = M 2 — 2MT + M-2k + 2 -M-1 + i)h - 

i=i 

M M fc 

—2 ^2 1 + 2 ^ 6i + 2 ^ 1 + 2(k + M — k)(X + T) 

i=l i=l i=l 

M 

= M(M + 1) - 2MT -2k- 2{M + 1) 6j - M(M + 1) + 

i=i 

M 

+2 ^ bi + 2A; + 2 MX + 2MT = -2{M + 1)X + 2X + 2 MX = 0, 
i=i 

and we obtain (m). □ 

Theorem 13. (Antisymmetry) 

For all M G N ; T G Z, ant/ rfeZ, 

seS(T,M+l,d) s'eS(M-T,M+l,-d) 

Proof. As in the proof of the preceeding proposition, we can match the states 
in the first sum with those in the second one. From property (Hi) in 12, we 
conclude that the classes, and thus the sum terms, are the same in each 
case. □ 
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Definition. For all d G Z, q = \¥ q \, M G N, and < T < M, let 
7 (d,T,M + l):= ^ ^(s) 

seS(T,M+l,d) 

be the asymptotic mass on all states with a given drain value d, at times 

EE (T,M+1). 

Theorem 14. For all d eZ, Me N, and < T < M, we /icwe 

7 (d, T, M + 1) = 7 (-d, M - T, M + 1). 
Proof. This follows immediately from Proposition 12 and Theorem 11. □ 
Definition. For < T < M and 1 < t < M + 1, let 

d(T,i) := /ioo(s) • 

seS(T,t) 

Also, let 

M 



d:= V d(T,M + l). 



M+ _ 

T=0 

Proposition 15. 

(i) For < T < M, d(T, M + 1) = -d(M — T, M + 1). 

(ii) For even M, we have d(M/2, M + 1) = 0. 

Proo/. (i) follows from Theorem 14, since d(T, M + 1) = J2 de %d->y(d, T,M + 
!))• 

(ii) follows from (i) with d(M/2, M + 1) — -d(M/2, M + 1) d(M/2, M + 
1) = 0. □ 

Theorem 16. For every M G N 7 d = 0. 

Proof. Using 15(i) (and 15 (ii) in case of even M), we have 

1 M _ 1 A M / 2 J_ 

5= ¥TT E^ M + 1 )= mTT E 3(T,M + l) + d(M-T,M + l) 

T=0 y T=0 

=0. □ 

Remark. In particular, Theorem 16 is an (aesthetical) reason to choose 
L a { n ) ~ r n ' a7+i"I ( an( i not L a (n) ~ n- j^j) as "typical" average behaviour, 
another reason is that for q — > 00 this same |~. . . ] value is the limit behaviour. 
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6. The Partition Model 



Definition. Let P M {K) G N, for M G N, K G N , be the number of parti- 
tions of K into at most M parts (equivalently, into parts of size at most M) . 

Definition. Let V(M,q) = Y,k=o P m{K) ' Q~ K - 

Proposition 17. a) The following initial values and recursion formulae hold: 
P X (K) = 1,VK G N, P M (1) = 1,P M (K) =0,\/K < 0,VM G N, and 
P M {K) = P M (K — M) + Pm-i(K). 
b) The generating function of Pm{K) in powers of g" 1 is 

OO M m 

q m 



V{M } q) = Y J PM(K)-q~ K =\{ 



K=0 



1?" 



R M-1 

c) Pm{K) & — — for fixed M and K — > oo. 

M!(M-1)! J J 

Proo/. See jH], Sections 2.5.10, 2.5.12 and 2.5.11. □ 

Remark. Observe that by c), for every K G No, we have only polynomially 
many states of class K, each with exponentially small probability q~ K • \i{sq). 
This leads to the concentration of mass on the states with small K. 

Definition. Let I m (s), 1 < m < M, count the number of actions / at battery 
m during a s (see Theorem 10). If K(s) = 0, put I m (s) = 0, 1 < m < M. 
Let (Ji, . . . , Im) be the ordered (Jj > Ii+i, 1 < i < M) version of {I m }. 

Corollary 18. Let jf^l be the number of inhibitions during all of the transi- 

M 

tions in a(s), similarly #iV<. Then I m = #1 = K(s). 

m=l 

Proof. The I m sum up to #J by definition. By Theorem 8, we have K(s) — 
K(s ) = #/-#jV<. With #A^< = and K(s ) = 0, #1 = K(s) follows. □ 

Theorem 19. 

(i) For 1 < M < 8, < T < M, 1 < t < M + I, and < K < 1200 - 100M, 
the state set S(T,t) contains exactly Pm(K) states with K(s) = K. 

(ii) For 1 < M < 8 and < K < 600 - 50M, fix a time (T ,t )- Then 
the . . . , Im) of all the Pm{K) states in 5 (To, to) with K{s) = K give the 
Pm{K) different partitions of K into M parts {including those of size 0). 

Proof. By numerical simulation over the mentioned ranges. □ 

Conjecture 20. The previous theorem holds for every M, T, t, K. 
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A graph for M — 2, showing states (61, 62, d) with K(s) < 4 and 
their associated partitions: 



Class K(s) 
T,t Part. [0,0] 



1 

[1,0] 



2 

[2,0] 



2 

[1,1] 



3 

[3,0] 



3 

[2,1] 



4 

[4,0] 



4 

[3,1] 



4 

[2,2] 



>4 

[*,*] T,t 



2,3 (-1,-1,0) (-1,0,-1) (0,-1,-1) (-2,0,0) (-2,-1,1) (0,-2,0) (-2,1,-1) (-1,-2,1) (0,0,-2) („) 2,3 




Accumulated Mass up to this column for 



q=2 

q=3 



Legend: 



0.375 
0.5925 



0.5625 
0.7901 



0.6562 
0.8559 



D I 



0.75 
0.9218 



N=,N< d-,b+ 



0.7968 
0.9437 



0.8437 0.8671 0.8906 0.9140 1.0 
0.9657 0.9730 0.9803 0.9876 1.0 
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Conjecture 21. (Theorem for M < 8) For every state s G S, 

q-K(s) 
= V(M,qY 

Proof. We assume the previous Theorem 19 or Conjecture 20. To normalize, 
we want to have 

1 = E c ° ■ *~ K{S) = c ° E Pa ^^) • ^ x = c ° p ^ M ' «)• 

With c := V(M,qY l = ^(sq), ^(s) := q~ K{s) /V(M, q) is a probabil- 
ity distribution (with E se5(T><) ^oo(s) = 1 for all < T < M, 1 < t < M+l), 
which is invariant under T. □ 



7. Asymptotic (n — » oo) Measure for the 
Linear Complexity Deviation 

Definition. Let the mass on states with drain (deviation) d be j(d,T,t) = 
^2ses{Ttd) ^oo(s), distinguished according to the timesteps (T,t). 

Numerical results indicate that 7 indeed depends only on the difference 
t-T: 

Theorem 22. For 1 < M < 8, < T < M, 1 < t < M + 1, and any finite 
field W q , let A := t — T . Then for every linear complexity deviation d G Z, 

7ld - p(m,,)^ ( » 9 <««)c-« ■n^'( ( *-DnS t *'-i: 

e_(A,/i) = /i(M-l + A)-g) 
£+ (A,/ i ) = h(A-h) + Q 
e (A,h) = min{e + (A,/i),e_(A,/0} 
depends only on the sign of d, and = means equality with precision at least 

^-(1200-100-M) 

Proof. By verifying all states with class up to 1200 — 100 ■ M in the partition 
model. The left and right side coincide up to precision g- 12 oo+iooA/ rj 

Remark. This involved about 2 39 or half a trillion states for M = 8. 
We used Victor Shoup's library NTL jH] (Thank you!). 
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Conjecture 23. For every MeN, < T < M, 1 < t < M + 1, and every 
finite field F g , with A and e(A, /z) as before, for every d G Z, we have exactly 

7 {d, t, t)= y ^ ^ — / } 

ti q (M+m + Mm [ii^ I+2 ( q k - i)j q h-(M+iy\ d \ 

the same formula as in Theorem 22, rearranged. 

Remark: The resulting values •y(d, T, M + 1) for M = 2 and M = 3, and 
d(T, M + 1) for M = 2, correspond with Niederreiter's and Wang's results in 
[TTl Thm. 3], [TTJ Thm. 4], and Thm. 11], resp., for n — > oo, see also 0. 
Observe that we use d = L — \n • |] , not L — n • §. 

8. The Law of the Logarithm 

We follow the approach by Niederreiter in [S] for the case M — 1. 

Theorem 24. For a// MeN, /or allO<T <M, and 1 < t < M + 1, t/iere 
exists a constant C(M,T,t) > (independent of d) such that: 

1 q -\dKM+i) < T ( djTjt ) < C(M,T,t) ■ q -\ d \< M + 1 \ 

1 (d,T,t) = e(q-^ M+ V). 



V(M,q) 

that is 



Proof. Lower bound: 

We distinguish cases d < 0, d > 0, and d = 0: 
a) < 

Let 6 := [^f^J and a :— —d — T — M • b that is 6 = -(a + d + T)/M. 
Then s* := (6, . . . , b, b + 1, . . . , b + 1, d; T, t) with b\ — ■ ■ ■ — bu-a = b and 
b M -a+i = ■ ■ ■ = b M = & + 1 is in S(T, t). 

The class of s* is #0*) = —n s + M ■ T + 2££ =1 (& + 1)(M + 1 - jfe) + 
2 J2k=a+i b ' \M + 1 — since after the sorting, the a batteries with value 
6+1 will be the largest, while d is the smallest value. 

Now, sorting starts with b, . . . , (d), . . . , b, b + 1, . . . , (d), . . . , 6 + 1, where d 
occupies the t-th place from the left, d moves to the right by M+ 1 — t moves, 
then all the 6's interchange with all the (6 + l)'s in a(M — a) transpositions, 
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yielding 7i> = M + 1 — t + a(M — a) and class 
K(s*) 

= - (M + 1 - t + aM - a 2 ) + MT + 26(M + l)M/2 

+ 2[(M + l)M/2 — (M — a + 1)(M — a)/2] 
= - M - 1 + t - aM + a 2 + MT - (d + T + a)(M + 1) 

+ (M + 1)M — (M — a + 1)(M — a) 
=|d|(M + l)-M-l+i-aM + a 2 + MT-MT-T-aM-a 

+ (M + 1)M - (M + 1)M - a 2 + 2aM + a 
= |rf|(M+l)-[(M-t) + l+T] 
<|d|(M + l) 

and already s* alone accounts for the lower bound. 

b) d > 

With 6, a, and s* as before, sorting now leads to a + a(M — a) transpo- 
sitions, since d goes to the left. As before, 

o+l M 

K(s*) = -a(M+l-a)+MT+2dM+2j2( b + 1 )( M + 1 - k )+ 2 H M + l ~ k ) 

k=2 k=a+2 

= -a(M+l-a)+MT+2dM+2b(M+l)M/2-2bM+2[M(M-l)/2-(M-a)(M-a-l)/2] 
= -a(M+l-a)+MT+2dM-(d+T+a)(M+l)+2(d+T+a)+M(M-l)-(M-a)(M-a-l) 
= -aM-a + a 2 + MT + 2dM - dM - MT — aM — d — T — a + 2d + 2T + 2a 

+M 2 - M - M 2 + 2a M - a 2 - M + a 
= dM + d + T + a - 2M = \d\(M + 1) - (M - T) - (M - a) < \d\(M + 1) 

c) d = 0: Lets* be the (only) state in S{T, t) with K(s*) = 0= |d|(M+l). 
Upper bound: 

We have 

7 (d,T,f)= ^ o- K(s) -Co 

seS(T,t,d) 

(use C := V(M,q)~ 1 , if you trust Conjecture 21), and 

M+l 

K(s) = -tt s + M • T + 2 • ^ 6 m • (M + 1 - m). 

m=l 
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a) d < 0: 

Let {bi, . . . , 6m} be ordered nonincreasingly as b[ > ■ ■ ■ > b' M (where ' 
indicates that d does not enter into the sort). 

We write the battery values as sum of their differences A k := b' M _ k — 
b' M _ k+1 > 0, f < k < M - 1 



M-m 



k=i 



(where for m = M the empty sum is zero). 

By the invariant d + T + J2 m b m = we must have 



M M-m 



m=l k=l 

M-1 



« M-b' M = -d-T- J2MM-k) (10) 

k=l 

When running A 1; . . . , A M _i through all values from N and setting b' M 
by (9-11) and then b' m by (8), we obtain all possible values for (b[, . . . ,b' M ) 
(and a lot more, since b M is taken from |dZ). 

Furthermore, each (6' 1; . . . , fo M ) corresponds to up to M! states (with dif- 
ferent permutation of the values) with {b' m } = {b m }. 

7i s can be bounded in general by < n s < ( M 2 ^ 1 ) (the maximum being 
attained in the case of b± < bi < • • • < b m ). 

The transition from (b[, ... ,b' M ),d to (b\, ... , &m+i) i-e. including d in 
the sort order, gives the inequality 

M M+l 

2 E ^ m ( M + l-m)+2-0-d<2^ b m {M + l-m) 

m=l m=l 

since at any rate a smaller b' m will be replaced by a larger or b' m _ 1 . 



22 



Putting things together, we have the upper bound 



MT x 



7 (d,T,f) < C ■■■ E M! ^ ( 2 K ^ 

AiGNq A M -i€No 



x ^2E: =1 -|(M + l-m)-2E: =1 
= _ g -2E^-l(M+l-™). C(M)T ) 

= - g -l d K M+1 ).C(M,T) 



where 

C(M,T) = C ]T ••• ]T Mix 

AieN Am-iGNq 



-MT-2T M - T -Sfli lA ^ M - fe ) 



is independent of d. 

b) In the case d > we follow the same idea, however we put d as first 
(largest) value of the sort order. We obtain 



m-2 



b' m = b' 2 -J2^k,2<k<M+l 



k=i 

with A k G No. 

The invariant d + T + J2 m b rn = requires 

M+l m-2 M-l 

m=2 fc=l fc=l 

_ r, = _ A _ T_ E£%(jM 
2 MM M 

Again up to M! states can be attached to one sorted tuple (& 2 , . . . , b' M+ i), 
again tt s < ( M 2 +1 ) , and introducing d (from the left) increases (if at all) the 
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values, i.e. b m >b' m . We obtain 



M+l 

-2^6 m (M+l-m) 

m=l 

M+l 

< -2-d- M -2^b' m {M +l-m) 

m=2 

= -2A« - 2 vf- A _ L + ESW _ £ A fe I (M + 1 - 



. m) 

. M M M 

m=2 \ k=l 
M+l , 

a 



-2Md + 2 V -— + C 
^ M 



m=2 



^2M + 2 ' M - 2 " tf ) + Cl 
-d(M+l) + Ci, 



where C% does not depend on d, and thus 

7 (d, T, t) < Co 53 ' ' • E M! ' ^ (Ma+1) ■ <T MT <r d(M+1)+Cl 

AiGNo A M _!GN 

= g -W(^+ 1 ).C(M,T). 

□ 

Lemma 25. (Borel-Cantelli) 

(i) Let Ai, A2, ■ ■ ■ be events which happen with probability ax, 0,2, ■ ■ ■ , resp. 
If now J^fceN ak < 00 > then with probability one only finitely many of the 

events Ak occur simultaneously. 

(ii) Let A\, A2, ... be independent events which happen with probability 
ax, 0,2, ■ ■ ■ , resp. 

If now X]fcGN a ^ = 00 > then with probability one infinitely many of the 
events Ak occur simultaneously. 

Proof. See Feller % VIII.3]. □ 
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Theorem 26. 

The Law of the Logarithm for Linear Complexity of Multisequences 

For all M e N and for almost all {in the sense of Haar measure on 
(F M ) ) sequences a G (F^) , we have 



lim sup 



and 



d a {n) 
logn 

lim inf < ^ a ^ n ^ -. 

n— >oo log Tl 



(M + l)logg 



(M+ l)logg' 



Proof We fix some e > and apply the Borel-Cantelli Lemma 31 (i) to the 
events 



d a (n) 



logn 



> 



1 + e 



(M+ l)logg' 



With L :-- 



log n 



(M+l) log g 



, the probability for A k is 



= ^ l(d, T,t)+J2 lid, T, t) < 2 • C(M, T) £ <T |d|(M+1) 



cZ=L 



= 2C(M,T)-_ 



-L(M+1) 



with accumulated probability 



^ 2C(M,T) x ~ . 

2^<73— (m+T)2^ 



i 11 (l + e) logn 
(M+l)- (M + 1)logq 



n=l 



n=l 



-(M+l)' 



2C(M,T) y> 

1 _ o-(M+l) Z_> 
n=l 



n- {1+£) < cx). 



For the inner bounds, we need independent events: 

Denote by rii, ri2, ■ ■ ■ the timesteps, when d = 0. If this sequence is finite, 
d — > — oo, since at least one battery no longer discharges. This event is of 
measure zero, requiring all discrepancies 5 pertaining to that battery equal 
to zero from some Uq on. 



Assume now an infinite sequence of these timesteps. Let L k :— 



log A: 



(M+l) log q 

and let A k be the event of (M + 1) • (L + 1) consecutive discrepancies, 
all zero, after n^. The events A k are independent with probability a& = 
q-(M+i)(L+i) ^ s j nce they belong to different, independent discrepancies. Now, 
within (L + 1)(M + 1) symbols, we have at least (L + columns and 
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thus [(L + — L actions oL (without intermediate discharges), 



M M+ 

and thus A)~ leads to d < —L. 
With 



n=l 



n=l 



n=l 



and Lemma 31 (ii), we get equality of the bounds. 



□ 



Corollary 27. For all MeN, for all e > 0, for almost all sequences from 
(Ff )°°, it holds 

y . r daiji) d a {n) 
—e < hm mi < hm sup < e 

n^oo n jwoo n 

Proof. Almost always, we have 



dJn) 



log n 



< 



1 



(M + l)logg 



d a (n) 



n 



< 



logn 



n(M + l)logg 



by the last theorem, and with 



1 



lim^ = 



(M + 1 ) log q n->oc n 
the statement follows. 
Theorem 28. With measure one, 



□ 



r ■ f L a {n) L a (n) 
hm mi = hm sup 



M 



n— >oo n 



n 



M + l 



Proof. From L a {n) = d a {n) + |"|^-] and the previous corollary, we have 



L a (n) d a (n) M 1, 
n n M+l ^n' 



and thus 



lim^ 



M 

hm h 



M 



n— >oo 77, n— >oo fj, 



M+l M+l 



□ 



In other words, we obtain again the result of Niederreiter and Wang [HIE] 



that 



Lg(n) 



M 
M+l 



with probability one, for all M G N. 
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9. Finite Strings 

Definition. For s E S, let the generation of state s be 



his) 



(M+l). 



M+l 

Conjecture 29. For every state s E S and every n E N 

0, n < g(s) A (T,t) ^ n 

A*n(s) = <( /ioo(s) • F(s), n = g(s) 

Hoo(s), n > g(s) A (T, t) = n 



with 



for 



M 



f(s)= n 



m=Mi(s) 



1' 



Mi(s) = M + 1 - # {1 < m < M | b rn = max{&i, ...,b M , d}} . 

In the case of the empty product for M 1 — M + l, F(s) = 1, and for 
M 1 (s) = l,F(s) = V(M,q). 

Conjecture 30. a) For g E N , let M) be the number of states that 
are reachable in the <?-th generation. Then 

b) #{sES\g(s)=g}=(°l? I )-(l I ). 



Conclusion 

We introduced the Battery-Discharge-Model BDM as a convenient container 
for all information about linear complexity deviations in (F^ 7 ) 00 . 

We obtained a closed formula for measures and averages for the linear 
complexity deviation, numerically proven for the cases M = 1,...,8, and 
conjectured for any M, which coincides with the results known before for 
M — 1,2, 3, but gives a better account of the inner structure of these mea- 
sures. In particular, the measure is a sum of M components of the form 

e(q-^ M+ V h ),h = l,...,M. 
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